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qq , We determine the three-neutrino mixing and mass parameters that are 

allowed by the solar and atmospheric neutrino data when vacuum oscillations 
are responsible for both phenomena. The global fit does not appreciably 
change the allowed regions for the parameters obtained from effective two- 
Q-i. neutrino fits. We discuss how measurements of the solar electron energy 

c~| , spectrum below 6.5 GeV in Super-Kamiokande and seasonal variations in the 

Super-Kamiokande, Ga, and BOREXINO experiments can distinguish the 
^ ■ different solar vacuum solutions. 
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I. INTRODUCTION 



Recent data from the Super-Kamiokande experiment PJ0] have strengthened the inter- 
pretation of the solar |3|-§| and atmospheric |||| neutrino anomalies in terms of neutrino 
oscillations. Oscillations have also been invoked to describe the appearance of electron neu- 
trinos and antineutrinos in the LSND experiment Because confirmation of the LSND 
results awaits future experiments and recent measurements in the KARMEN detector ex- 



clude part of the LSND allowed region ||11||, a conservative approach is to assume that 



oscillations need only account for the solar and atmospheric data. Then the two mass- 
squared difference scales in a three-neutrino model are sufficient to describe the data. A 
number of three- neutrino models have been investigated |12|- |17|1 . One attractive possibility 
is that both the atmospheric and solar v e oscillate maximally [0] or near-maximally [T5| 
at the Sm 2 tm and 5m 2 un scales, respectively. 

Most detailed oscillation fits have been done separately for the solar and atmospheric data 



in effective two-neutrino approximations [18|. In this paper we make full three-neutrino fits 



to the solar and atmospheric oscillation data to determine the allowed values for the general 
three-neutrino mixing matrix under the assumption that one mass-squared difference, 5m r 
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explains the atmospheric neutrino oscillations and that the other independent mass-squared 
difference, 5m 2 un <C 5m 2 tm , explains the solar neutrino oscillations via the vacuum long- 



wavelength scenario |19| , p0[ . We choose a particular parametrization for the three-neutrino 
mixing in which the expressions for the atmospheric and solar neutrino oscillations share 
only a single parameter, and discuss the degree to which the two phenomena decouple. We 
find that the extension from two to three neutrino species does not improve the separate fits 
to either the atmospheric or solar data. Although pure — > v T oscillations of atmospheric 
neutrinos are favored, there exist three-neutrino solutions with non-negligible v„ <-> v e os- 



cillations, even with the constraints from the CHOOZ reactor experiment ||21|| . Hence both 
Vfj, <-> v e and v e — > v T oscillations may be observable in future long-baseline experiments. 
We also find that future measurements of seasonal variations in the solar neutrino signal can 
provide an unmistakable sign of vacuum oscillations, and can further constrain the allowed 
parameter regions. 

In Sec. II we review the formalism for oscillations of three neutrinos relevant to the 
atmospheric and solar phenomena. In Sec. Ill we evaluate the current allowed two-neutrino 
parameter regions, and briefly review the evidence indicating that two distinct 5m 2 are 
needed to describe the solar and atmospheric data. In Sec. IV we obtain three-neutrino 
solutions from a combined fit to the atmospheric and solar data. We conclude with a 
discussion of future tests of vacuum oscillations in Sec. V. 

II. OSCILLATION ANALYSIS 
A. General probability expressions 

The survival probability for a given neutrino flavor v a in a vacuum is 

P{v a -> v a ) = 1 -4£ \U a /\U ak \ 2 sm 2 A ife , (1) 
k<j 
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where U is the neutrino mixing matrix (in the basis where the charged-lepton mass matrix is 
diagonal), A jk = 5m 2 k L/AE = 1. 27 (5m 2 k /eV 2 )(L /km) /(E/GeV), 5m 2 jk = m 2 -m 2 k} and the 
sum is over all j and k, subject to k < j. The matrix elements U a j are the mixings between 
the flavor (a = e, //, r) and the mass (j = 1, 2, 3) eigenstates, and we assume without loss of 
generality that mi < < m-&. The solar oscillations are driven by IA21I = A sun and the 
atmospheric oscillations are driven by |A 3 i| ~ |A 32 | = A atm 3> A t 



The off-diagonal vacuum oscillation probabilities in a three-neutrino model are 

P(u e -»• = 4 |f/ e3 f/; 3 r sin 2 A atm - 4Re{U el U^U^U^} sin 2 A sun - 2 Jsin 2A sun , (2) 
P(u e ^u T )=A \U e3 U* rZ \ 2 sin 2 A atm - ARe{U el U* e2 U* x U T2 } sin 2 A sun + 2 Jsin2A sun , (3) 
P(^ -»• i/r) = 4 |^; 3 | 2 sin 2 A atm - ARe{U,aU; 2 U* Tl U T2 } sin 2 A sun - 2 J sin 2A sun , (4) 

where the CP-violating "Jarlskog invariant" p4[ is J = J2k,-y e ijk^ai3-yInT'{U a iU*jUp i Ui3j} for 
any a, j3, i, and j (e.g., J = Im{U e 2U* 3 U* 2 U^3} for a — e, (5 — fi, i — 2, and j = 3). The CP- 
odd term changes sign under reversal of the oscillating flavors, or if neutrinos are replaced 
by anti-neutrinos. We note that the CP- violating probability at the atmospheric scale is 
suppressed to order 5m 2 un / '8m 2 atm1 with the leading term cancelling in the sum over the two 
light-mass states; thus, P(u a — > Vp) = P{p a — > Vp) (and therefore P{y a — > vp) = P{yp — > v a ) 
from CPT invariance) at the atmospheric scale. 

Without loss of generality we work in the basis where the charged-lepton mass matrix 
is diagonal. The matrix U that relates the flavor eigenstates to the mass eigenstates may 
be parametrized in terms of three Euler angles and one (three) phases for Dirac (Majorana) 
neutrinos. This can be understood as follows. In the Dirac case, U is analogous to the CKM 
mixing matrix in the quark sector, in which there are three Euler angles and one phase. For 
Majorana neutrinos, the 3x3 mass matrix in the flavor basis must be symmetric but may 



be complex [ 55f] . There are 12 independent parameters, which can be taken as 3 real mass 
eigenvalues, 3 real Euler angles, and 6 phases. Then three of the phases can be absorbed into 
the definitions of the fields. However, only one of the three remaining phases has physical 
consequences in neutrino oscillations. Therefore for either Dirac or Majorana neutrinos we 
can choose the following parametrization for U [EH] 



(5) 
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sm9j and 8 is the CP-violating phase. 



B. Atmospheric and long-baseline experiments 

The experimental indications are that 5m 2 atm ~ 10 -3 eV 2 and that Sm 2 un ~ 10~ 10 eV 2 
for a vacuum oscillation explanation of the solar neutrino data. Then for the oscillations of 
neutrinos in atmospheric and long-baseline experiments with L/E> 10 2 km/GeV, the A sun 
terms are negligible and the relevant oscillation probabilities are 

PK *V) = 1 - (4 sin2 2 °2 + sj sin 2 26 x ) sin 2 A atm , (6) 
P(u e v e ) = 1 - sin 2 29 1 sin 2 A atm . (7) 
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P(p„ 



s\ sin 2 26»i sin 2 A atm . 
cl sin 2 26»i sin 2 A atm , 
c\ sin 2 2# 2 sin 2 A 



atm 



(8) 
(9) 
(10) 

When 6*i = (i.e., U e ^, = 0), these reduce to pure — > u T oscillations with amplitude 
sin 2 2^2, and v e does not oscillate in atmospheric and long-baseline experiments. 

We define the oscillation amplitudes A^f m , A e a i m i Aatm: A e J tm , and A^J tm as the coefficients 
of the sin 2 A atm terms in Eqs. (p|)-(|Tt]|), respectively. The neutrino parameters can then be 
determined from the atmospheric neutrino data by the relations 



NJN° = a [(1 - (S) A±) + r (S) A% m 



(11) 



and 



a 



(1 — {&) Aatm) + r (S) A^trn 



(12) 



where N° and N° are the expected numbers of atmospheric e and fi events, respectively, 
r = N°/N°, (S) is sin 2 A atm appropriately averaged, and a is an overall neutrino flux 
normalization, which we allow to vary following the SuperK analysis 0. 

SuperK presented N^/N" and N e /N° for eight different L/E bins from a 535 day 
exposure. The data were obtained by inferring an L/E value for each event from the zenith 
angle 9e and energy of the observed charged lepton Eg and comparing it to expectations 



from a monte carlo simulation based on the atmospheric neutrino spectrum |27|] folded with 
the differential cross section. Due to the fact that the charged lepton energy and direction 
in general differ from the corresponding values for the incident neutrino (or antineutrino), 
the L/E distribution involves substantial smearing. We estimate this smearing by a monte 



carlo integration over the neutrino angle and energy spectrum ]28| weighted by the deep- 
inelastic differential cross section [2~9|. We generate events with E v and 6 V , and determine 
the corresponding Ei and 9e for the charged lepton. We bin the events in L/E u , using $£ 
to determine L and an estimated neutrino energy inferred from the average ratio of lepton 
momentum to neutrino energy, E^ = Eg {E u /Eg), analogous to the SuperK analysis 0. We 
calculate a value for (sin 2 A atm ) for each L/E bin for a given value of 5m 2 atm . We can then 



fit Eqs. ( [TT]) and (|I2"D to the data and determine the neutrino parameters in Eqs. (||)-(§D 
and the normalization a. Without loss of generality we take Sm 2 atm to be positive. 

favored by our fits, 



We do not consider matter effects in our analysis. For the 5m 2 atrn 
matter effects are small for the sub-GeV neutrinos that constitute most of the data 
Also, as evidenced by our fits, the dominant oscillation is — > u T , which is not greatly 

However, matter effects could be important for the neutrino flux 



affected by matter |31| 
with smaller 5m 2 atm /E : 



i.e., multi-GeV data in solutions with 5m 2 



atm 



< 



10 3 eV , or in 



long-baseline experiments 3C . 



C. Solar experiments 



For neutrinos from the sun L/E ~ 10 km/GeV, and the sin A atm terms oscillate very 
rapidly, averaging to h. Then the oscillation probabilities are 
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P(u„ 



P(v P 



P{v P 



v e ) = 1 - - sin 2 26>i - c\ sin 2 26> 3 sin 2 A sun , 

^m) = sin2 2 #i + 4c 2 s 3 c 3 s 3 c 3 (c 2 - s 2 s 2 ) + sis 2 c 2 cos 26> 3 cos 5 
-2sic 2 s 2 c 2 s 3 c 3 sin 5 sin 2A swl , 



sin 2 A, 



u T ) = -c\ sin 2 20i + 4c 2 s 3 c 3 



S 3 C 3( S 2 — s \ c \) ~ s i s 2C2 cos 26* 3 cos 5 sin 2 A 



+2sic 2 S2C2S 3 c 3 sin 5 sin 2 A s 



(13) 



(14) 



(15) 



where P(v a — > z^g) = P (z?g — > P Q ) from CPT invariance and P{p@ — > f a ) m &y be found 
from P[y a —>■ vp) by changing the sign of 5. Since only the sum of oscillation channels 
P[y e — > Vy) + P(^ e — ► ^r) = 1 ~~ P{y e — ► i/ e ) is tested in solar experiments, the CP-violating 
parameter J cannot be constrained from solar measurements. To see the effects of CP 
violation one must measure 



P(v a -> i/ e ) - P(u e -> = P(z/ e -> i/ r ) - P(z/ T -> 

= P(z/ r -> Vp) - P(z/ M -> i/ T 
= 2siC 2 s 2 c 2 s 3 c 3 sin 5 sin 2A S 



(16) 
(17) 
(18) 



the corresponding differences for antineutrinos (which have the opposite sign), or combina- 
tions that violate CP explicitly, such as P(v n — > f e ) — Pipu —* v e ). 

When 61 = (i.e., ?7 e3 = 0), the solar oscillation acts like a simple two- neutrino oscillation 
with amplitude sin 2 2# 3 ; although the oscillations of v e may involve both and u T , the 
individual channels u e — > z/ M and z/ e — > i/ r are not measurable in solar experiments, since 
the neutrino energies are below the thresholds for fx and r production. The parameter Q\ 
measures the extent to which the solar oscillations have a constant component coming from 
the atmospheric oscillation scale [[32[ . 

In our solar fits we use the shape of the standard solar model (SSM) neutrino fluxes and 
absorption cross sections in Ref. |6j and the new normalizations in Ref. j7|. For the 37 C1 and 
71 Ga cases, we fold the neutrino oscillation probability with the neutrino absorption cross 
section and expected neutrino flux to obtain the expected number of events, which can then 
be compared to the expectations of the standard solar model. We also allow an arbitrary 
normalization factor /3 for the 8 B neutrino flux. The expected number of neutrino events is 
then 



N = J CiP(lS e -> V e )((3(f)B + <pnon-B)dE u , 



(19) 



For the Super-Kamiokande case, for which the interaction is ve — > ve, the events are 
binned by outgoing electron energy, and we have included the effects of the detector resolu- 
tion. The number of events per unit of electron energy is 



dN 




P(y 9 ^ Vo ) + ^[l-P( Ve 



,)\\G{E' e ,E e )4> B dE v dE' e (20) 



where dacc/dE e (da^c I dE e ) is the charge-current (neutral-current) differential cross section 
for an incident neutrino of energy E v and G(E' e , E e ) is the probability that an electron of 
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energy E' e is measured as having energy E e ; the electron energy resolution is taken from 
Ref. 0. The events are then put into 0.5 GeV bins starting at the 6.5 GeV threshhold for 
the detector, to compare with the Super-Kamiokande data. We take as the input solar data 
the Homestake 37 C1 rate [g, the GALLEX and SAGE f| 71 Ga rates, and 16 bins of the 
Super-Kamiokande detected electron energy spectrum We then fit Eqs. ( [19| ) and ( |20|) 
to the data and determine the neutrino parameters in Eq. (|13|) and the 8 B neutrino flux 
normalization j3. 

In our solar fits we combine day and night results, since they should not be appreciably 
different for vacuum long-wavelength oscillations. We compare the model predictions with 
the time-averaged data, since with current statistics the Super-Kamiokande data does not 
reflect any seasonal variation [[l]. The possibility of detecting a seasonal variation in Super- 
Kamiokande and other experiments with improved statistics is discussed in Sec. V. 



III. TWO-NEUTRINO SOLUTIONS 
A. Independent solutions for atmospheric and solar data 



In the limit 9\ = (i.e., U e ^, = 0), the atmospheric and solar oscillations decouple p2|j33 



and effectively reduce to two separate two-neutrino solutions, each with its own oscillation 
amplitude and mass-squared difference. Also, the only oscillation channel for atmospheric 
neutrinos is —>■ v T . More specifically, for atmospheric neutrinos 

PW ~* v v) = 1 " sin2 2 ^2 sin 2 A atm , (21) 
P(u e -> u e ) = 1 , (22) 

P(v e "v ft )=0, (23) 

and for solar neutrinos 

P{y e -> u e ) = 1 - sin 2 28 3 sin 2 A sun . (24) 
In the U e s = limit, fits to the atmospheric and solar data may be made independently. 



B. Atmospheric data 

For the two-neutrino case with 9\ = (only — ► v T oscillations) our best fit parameters 



are 



5m 2 atm = 2.8 x 10- 



sin 2 26 2 
a 



eV 2 



1.00 
1.16 



(25) 
(26) 
(27) 



where a is the overall flux normalization in Eqs. (11) and (12), with Xmin = 7.1 for 13 degrees 
of freedom (DOF). This \ 2 /DOF corresponds to a goodness-of-fit of 90%. The 95% C.L. 
allowed region for 5m 2 atm versus sin 2 2# 2 is shown in Fig. |I[ Our result is very similar to the 
fit obtained by the SuperK collaboration. If we set a = 1 (i.e., assume that the theoretical 
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flux normalization is exact), the best fit has x 2 /DOF = 22.8/14, which is acceptable only 



at the 6% CL. The calculated flux has a normalization uncertainty of about ±20% [34 . 
As reported by the SuperK collaboration, the other two-neutrino case with pure vacuum 
<-> v e oscillations (which corresponds to 9 2 = vr/2) does not give a good fit to the data 

(possible effects of matter are discussed at the end of Sec. II. B). We find that the minimum 

X 2 /DOF for this case is 81.9/13, corresponding to a goodness-of-fit of 4 x 10 -12 . Therefore 

the <r-¥ u e scenario for atmospheric neutrinos is strongly disfavored by the SuperK data. 

The results of the two-neutrino fits to the atmospheric data are summarized in Table |. 

Large amplitude — > v e oscillations are also excluded by the CHOOZ reactor data for 

Sm 2 atm > 10^ eV 2 . 



C. Solar data 

For the effective two- neutrino oscillation formula with d\ = 0, our best fit to the combined 
solar data yields the parameters 

5m 2 sun = 7.5 x 10- 11 eV 2 , (28) 
sin 2 2^ = 0.91, (29) 
0= 1.62, (30) 

with XmiJ DOF = 21.6/16, acceptable at the 16% C.L. The 95% C.L. allowed regions for 
Sm 2 tm versus sin 2 2# 3 are shown in Fig. We note that our allowed regions are very similar 
to those obtained in Ref. |$5| with a somewhat different analysis. Taken at face value, the 
preferred value of the 8 B normalization would suggest that the SSM underestimates the 8 B 
neutrino flux by a sizable amount. However, this seems unlikely because most alternative 
solar models give a lower 8 B neutrino flux. The best fit values for (3 = 1 (the SSM 8 B 
spectrum normalization) are 

5m 2 uri = 6.5xlO- n eV 2 , (31) 
sin 2 29 3 = 0.74 , (32) 

with XmiJDOF = 26.5/17, acceptable at the 7% C.L. 

We have also searched for the best-fit to the solar data in each of the 5m 2 un "finger" 
regions in Fig. |2|. The results are given in Table O, where these four best fits have been 
labeled A, B, C, D in order of ascending 5m 2 un . These fits correspond to vacuum oscillation 
wavelengths (for a typical 8 B neutrino energy) of approximately |D, |D, |Z), and \D, 
where D is the Earth-Sun distance. We see that the best fits in the higher 5m 2 sun regions 
all have a lower 8 B normalization. They also tend to fit the Super-K data better, but do 
somewhat worse on the radiochemical experiments. Interestingly, the two solutions with the 
largest 5m 2 un have sin 2 2# 3 ~ I. 

Because the apparent suppression in the 37 C1 measurement differs from that of the 71 Ga 
and SuperK cases, we consider the possibility that it had an unknown systematic error. The 
results of fitting to just the 71 Ga and SuperK data are listed in the lower half of Table 0. 
The best fit now occurs for a larger 5m 2 un , and the goodness-of-fit improves. The allowed 
regions do not qualitatively change. However, for the larger Sm 2 un regions, the 71 Ga data is 
more easily accomodated by small changes of the parameters from their values in the global 
fit. 
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D. Existence of separate mass scales 



In making our three-neutrino fits we assume that separate mass-squared difference scales 
are necessary to explain the atmospheric and solar neutrino data. If the 5m 2 scales were not 
distinct, or if one of the mass-squared differences were used to explain the LSND data, then 
either the solar or atmospheric probabilities would be in a region of L/E where the oscil- 
lations have averaged, and there would be no energy dependence. An energy-independent 
suppression due to oscillations would be equivalent to letting the overall normalization vary. 
Oscillation scenarios where there is not a separate mass scale associated with solar neutrinos 
have been considered |3~61 . 



It has already been demonstrated in the literature that an energy-independent suppres- 
sion is strongly disfavored by the solar data |37|. Even if one ignores the 37 C1 data, and 



assumes that the 71 Ga and ve experimental rates are both consistent with an overall sup- 
pression by 50%, the spectrum distortion of 8 B neutrinos measured in the ve experiments 
disfavors an energy-independent suppression. We have updated this analysis, allowing for 
an overall flux suppression in addition to the variation of the 8 B neutrino normalization. 
The lowest value P(v e — > v e ) can achieve with three neutrinos when all oscillations are av- 
eraged is |. We found the best fit to the solar data with overall depletion between | and 
unity and arbitrary 8 B normalization has ^/DOF = 48.1/17, which is ruled out at the 
99.99% C.L. If the 37 C1 data are ignored, the best fit has only x 2 /DOF = 25.0/16, ruled 
out at the 93% C.L. Therefore the distortion of the electron energy spectrum present in the 
solar neutrino data favors the existence of a separate mass scale for the oscillation of the 
solar neutrinos, which justifies the form of Eqs. (|^)-([T0|) and (|i~3|)-(|l5|). 

If one 5m 2 is used to describe the LSND data and the other is used to describe the solar 
data, then the oscillation due to 8m 2 LSND will be averaged for the L/E of the atmospheric 
neutrinos, so that the oscillation probabilities are independent of both energy and zenith 
angle. We find the best fit in this scenario to have x 2 /DOF = 33.2/13, which is excluded 
at the 99.8% C.L. 



IV. THREE-NEUTRINO SOLUTIONS 



A. Bi-maximal solution 



The atmospheric data favor maximal mixing of atmospheric with v T and no mixing 
with v e . The solar data also suggest, although not as strongly, that solar neutrinos may also 
mix maximally, or nearly maximally. If we require both atmospheric and solar oscillations 
to be maximal, there is a unique three-neutrino solution to the neutrino mixing matrix ||14|| , 
which corresponds to Q\ = and 6 2 
for atmospheric neutrinos are 



#3 = 7r/4. The corresponding oscillation probabilities 



P(v P 



"HI 
V P ) 



1 
1 





sin 2 A 



atm i 



(33) 
(34) 
(35) 



and for solar neutrinos 



S 



P[y e -> v e ) = 1 - sin 2 A sun . , (36) 
P[y e -> = P(u e -> v r ) = - sin 2 A sun . (37) 

One interesting aspect of this solution is that the solar u e oscillations are 50% into and 
50% into z/ r , although the flavor content of the v e oscillation is not observable in solar 
experiments. Further properties of the bi-maximal and nearly bi-maximal solutions are 



discussed in Ref. 14 . 



B. Atmospheric data 

A full three-neutrino fit to the atmospheric data with one 5m 2 atm scale has been made in 
Refs. [23,38]. In terms of the oscillation parameters defined in Sec. II, our best fit values for 



the four parameters are 

5m 2 atm = 2.8 x KT 3 eV 2 , (38) 

sin 01 = 0.00, (39) 

sin 2 20 2 = 1.00, (40) 

a = 1.16, (41) 

with X 2 miJ D OF = 7.1/12, acceptable at the 85% C.L. In Fig. | we show the 95% C.L. 
allowed region for sin 2 2# 2 versus Sm 2 tm for various values of sin 6i when the flux normaliza- 
tion a is allowed to vary. Although sin #i = is favored, nonzero values are allowed, which 
permit some <-> v e and v e — > v T oscillations of atmospheric neutrinos. In Fig. f| we show 
the 95% C.L. allowed region for sin#! versus 8m 2 tm when a and sin 2 2# 2 are allowed to vary. 



Another limit on sin#! comes from the CHOOZ reactor experiment |£TJ that measures 
v e disappearance 

Alt = 4P e3 (l - P e3 ) £0.2, (42) 

which applies for 5m 2 atm ^ 2 x 10 -3 eV 2 . The exact limit on varies with 5m 2 tm , and 
for 5m 2 atm < 10~ 3 eV 2 there is no limit at all. For 5m 2 atm = 2.8 x 10~ 3 eV 2 and a = 1.16, 
sin 6*i is constrained to be less than 0.19. The result of imposing the CHOOZ constraint is 
shown in Figs. |3| and [|. In Fig. |4] we see that the range of sin 9\ allowed by the fit to the 
atmospheric neutrino data and the CHOOZ constraint is 

< sin 0i < 0.29, (43) 

at 95% C.L. 



C. Solar data 

A full three-neutrino fit to solar data can be made by varying 9±, #3, Sm 2 un and the 8 B 



flux normalization f3, and using the expression in Eq. (13) for the oscillation probabilities in 



Eqs. ( |T9| ) and (p0|). We find the best fit values for the solar parameters 
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Sm 2 atm = 7.5 x 10- 11 eV 2 , (44) 

sin0i = O.OO, (45) 

sin 2 26 3 = 0.91 , (46) 

= 1.62, (47) 

with ■)^ nin /DOF = 21.6/15, acceptable at the 12% C.L. If the 8 B normalization is fixed at 
the SSM value (f3 = 1), the best fit is 

5m 2 tm = 6.5xlO- n eV 2 , (48) 
sin0i = O.OO, (49) 
sin 2 2^3 = 0.74, (50) 

with xLin/ DOF = 26.5/16, acceptable at the 5% C.L. The addition of the extra parameter 
sin#! does not improve the fit to the solar data. In Fig. |^ we show the 95% C.L. allowed 
region for sin 2 2# 3 versus Sm 2 un for various values of sin 9\ when (3 is allowed to vary. As in 
the atmospheric neutrino case, sin 6\ = is favored, although non-negligible nonzero values 
are allowed. In Fig. [| we show the 95% C.L. allowed region for sin#i versus Sm 2 un when f3 
and sin 2 26^ are allowed to vary. The range of sin 6\ allowed by the solar data is 

< sin < 0.49, (51) 

at 95% C.L. 



D. Global atmospheric and solar fit 

Since the atmospheric and solar fits have only one common parameter, sin^, and the 
best separate fits to the two data sets both have sin^i = 0, then the best fit to the combined 



solar and atmospheric data sets can be identified immediately as given by Eqs. (|3"8|)-(f4~T|) 
and (H)-©, with XmijDOF = 28.7/28, acceptable at the 43% C.L. To see how the fits 
vary with the common parameter, in Fig. [7| we show Xmin f° r ^he s °l a r and atmospheric data 
sets versus sin#i, as well as the combined x 2 ■ From the figure we see that the dependence 
of Xmin on sin^x is relatively weak for sin 6^ < 0.3 for the solar data set. Therefore, the 
parameter regions allowed by the combined solar and atmospheric data sets are essentially 
the most stringent of those allowed by the solar and atmospheric data sets separately. 



V. DISCUSSION 



A. Long-baseline oscillations 



The MINOS [39], K2K [40], ICARUS [41] and NOE [41] experiments can test for v. 



oscillations for 5m 2 atm > 10 



eV , and MINOS and ICARUS can also search for v u 



atm " 



Together these measurements could more precisely determine sin 2#2, sin#i, and 5m\ 
Further measurements of atmospheric neutrinos will also help constrain these parameters. 
Full three-neutrino fits including the solar neutrino data |32| can then determine one of the 
remaining two independent parameters in the mixing matrix, e.g., sin 2 2# 3 , and Sm 2 



,2 

sun' 
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There is new physics predicted when sin#i 7^ 0, i.e., v e — > z/ r oscillations with leading 
probability given by Eq. (||). The allowed range at 95% C.L. for the v e — > v T oscillation 
amplitude A e J tm versus 8m 2 atm is depicted in Fig. || the effect of the CHOOZ constraint is also 



shown. The maximal v e — > v T amplitude of about 0.15 occurs at 8m atm = 1.7 x 10 eV . 
These v e —>■ v T oscillations could be observed by long-baseline neutrino experiments with 
proposed high intensity muon sources |f42|-fi4|1 , which can also make precise measurements 
of z/ M <-> v e and z/ M — > v T oscillations. Sensitivity to A e J tm (5m 2 atm / 'eV 2 ) 2 > 2.5 x 10~ 9 is 
expected j|2] for the parameter ranges of interest here; matter effects would not be large 



for an experiment from Fermilab to Soudan pD[ . For 5m 2 atm = 2.8 x 10~ 3 eV 2 , A e J tm could 
be measured down to 3 x 10~ 4 ; the expected sensitivity versus 5m 2 atm is shown in Fig. |8|. 
Precise measurement of v e — > v T and «-> v e oscillations in such a long-baseline experiment 
would uniquely specify the CP-conserving part of the three-neutrino model. 



B. Future solar tests 

In this section we discuss how future solar neutrino measurements can distinguish be- 
tween the different solar vacuum oscillation solutions. We consider only the two-neutrino 
solutions since the discussion easily generalizes to the three-neutrino case. 

One interesting feature of vacuum long-wavelength solutions is that they can give a rise 
in the fraction of surviving z/ e 's for higher electron energies, in agreement with the SuperK 
measurement |IJ. Figure [] shows the ratio of the electron energy spectrum to the SSM 
prediction for two different vacuum long-wavelength scenarios and the SuperK data. Future 
measurements at SuperK will improve the statistics in the high E e bins. However, with an 
hep flux that is ~ 25 times greater than the SSM result, this rise at high E e could be due 
to contributions of hep neutrinos |[45|| . SuperK is also planning to lower their threshhold for 
electron detection to 5 MeV; this is particularly important since the number of events rises 



for energies just below the current threshhold of 6.5 MeV. In Fig. [10] we show predictions for 
the electron energy spectrum in SuperK in the range 5 MeV < E e < 10 MeV for the four 
oscillation solutions A, B, C, and D in Table [0]. It is evident that Solutions A and B may 
be distinguishable from C and D, and from each other, using data at lower E e . The energy 



spectrum measured in the SNO experiment |46| will also provide additional information on 



the energy dependence of the spectrum suppression. 

Another feature of vacuum solar solutions is that they may cause a detectable seasonal 
variation as the distance between the Earth and Sun varies P?| , fi8"|] . We define two seasonal 
asymmetry parameters, 

Ax = W"-*s) (52) 

(N w + N S + N F + N SP ) ' V ; 

= N W + N S -N F - N SP 
2 Nw + Ns + Nf + Nsp 1 [ ' 

where Nw, Ngp, N$, and Np are the number of events collected in the time periods from 
November 20 to February 19 (winter, Earth closest to the Sun), February 20 to May 21 
(spring), May 22 to August 20 (summer, Earth farthest from the Sun), and August 21 to 
November 19 (fall), respectively. Since N F = N SP (the same range of distances is covered), 



11 



A\ and A2 are the only independent quantities that may be constructed from the four 
seasonal measurements. The quantity Ai is similar to the seasonal asymmetry defined in 
Ref. |[49||, and A\ and A 2 are similar to the first two harmonics in the analysis of the first 



paper in Ref. J48[. The parameter A\ has a significant nonzero value when the oscillation 
probability increases or decreases monotonically as the Earth moves from perihelion to 
aphelion, while A2 has a significant nonzero value when the oscillation probability reaches 
a local extremum somewhere between perihelion and aphelion. Due to the 1/r 2 dependence 
of the solar neutrino flux and the 3.3% change in r from perihelion to aphelion, A\ and A 2 
have the values 0.030 and 0, respectively, in the absence of oscillations. The asymmetry A\ 
is strongly correlated with the electron energy spectrum distortion; such a correlation is a 



distinctive characteristic of vacuum oscillations 49 



For 8 B neutrinos, due to their relatively higher energies and longer oscillation wave- 
lengths, the seasonal variation is less pronounced than for solar neutrinos of lower energy. 
The total SuperK event rate as a fraction of the SSM value versus time of year |]J is shown 
in Fig. [11]; also shown is the prediction for Solution A in Table and the prediction for 
no oscillations. The results for other oscillation solutions are similar to the prediction for 
Solution A. Oscillation solutions produce an enhanced seasonal effect |47|-449|. The SuperK 
experiment has not observed a significant seasonal variation with the current data sam- 
ple, but could be sensitive to the effects predicted by vacuum oscillations with increased 
statistics. 

In order to extract the most information from the SuperK data it is advantageous to 
plot the seasonal asymmetries versus the observed electron energy. In Fig. [12] we show A\ 
versus E e for Solutions A, B, C, and D of Table P]. Although the asymmetries are not large, 
each solution has a characteristic shape, especially Solution A. The energy dependence of 
Ai clearly distinguishes the oscillation scenarios from the asymmetry induced only by the 
seasonal flux variation. Similar measurements can be done in the SNO experiment [f46 



The deviation of A2 from the value for no oscillations is at most about 0.003, and does not 
provide a good discrimination between models. 

The GALLEX and SAGE 71 Ga experiments may also exhibit a seasonal variation with 
increased statistics if there are vacuum oscillations of solar neutrinos pD , pT| . In Fig. [13] we 
plot the predictions for A 2 versus Ai in the GALLEX and SAGE experiments for a range 
of 5m? sun from each of the four "finger" regions in Fig. |2|. Most of the seasonal asymmetry 
in the 71 Ga experiments is due to the monoenergetic 7 Be neutrinos which constitute about 
25% of the signal in the SSM. 

The BOREXINO experiment [52] will primarily measure the 7 Be neutrinos (E u = 0.862 
MeV) using the process ve — > v e . The final-state electron kinetic energy T e has a maximum 
value of 0.665 MeV for 7 Be neutrinos. The pp neutrinos have a maximum T e of 0.26 MeV, 
and also there is a considerable background for T e < 0.25 MeV. Therefore selecting events 
in the range 0.26 MeV < T e < 0.665 MeV will preferentially select the 7 Be neutrino signal. 
Neutrinos from the pep reaction and the CNO cycle will also give final-state electrons in this 
energy range, but the 7 Be neutrinos represent more than 80% of the signal, assuming the 
SSM. 

Predictions for event rate and seasonal variation in the BOREXINO experiment for var- 



ious vacuum oscillation parameters have previously been discussed in the literature |[47| , |50 



Here we examine how the seasonal asymmetries in Eqs. (|52|) and (^) may be used to fur- 
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ther discriminate the different vacuum neutrino solutions in Table ||. In Fig. [14] we plot the 
predictions for A 2 versus Ai in the BOREXINO experiment for a range of bm 2 sun from each 
of the four "finger" regions in Fig. 0. The seasonal asymmetries are potentially larger than 
in the 71 Ga case since the 7 Be neutrinos are a much larger fraction of the signal. Although 
there are some regions where the predictions of two or more solutions overlap, combining the 
asymmetry information with the event rate should significantly reduce the allowed parameter 
regions for the vacuum solutions, and could select the appropriate "finger" region. 
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TABLE I. Best-fit two-neutrino solutions to the atmospheric data for different oscillation sce- 
narios. The sin 2 29 in each case corresponds to the effective two-neutrino oscillation amplitude, 
with 9 = 82 (81) for v„ — > v T (v„ — > v e ) oscillations. 



Oscillation 


5m 2 atm (10" 3 eV 2 ) 


sin 2 29 


a 


XL/DOF 


Goodness-of-fit 


Vfi -> I/ T 


2.8 


1.00 


1.16 


7.1/13 


90% 


^ -> 


1.1 


0.85 


1.00 (fixed) 


22.8/14 


6% 


^ -> f e 


1.1 


0.88 


0.71 


81.9/13 


4 x 10- 12 



TABLE II. Best-fit two-neutrino solutions to the solar data in the four "finger" regions of 
<5m 2 un depicted in Fig. ^. The top half of the table corresponds to a fit to all solar data, while the 
bottom half shows the results for a fit with the 37 C1 data excluded. The x 2 sums do not add to 
the total in some cases due to rounding. 



Solution 


5ml un (10- 10 eV 2 ) 


sin 2 26» 3 


P 


37 C1 


71 Ga 


Super-K 


Xtot 


Goodness-of-fit 


A 


0.75 


0.91 


1.62 


0.8 


1.5 


19.4 


21.6 


16% 


B 


2.49 


0.86 


0.84 


5.9 


2.2 


18.2 


26.3 


5% 


C 


4.40 


0.97 


0.80 


6.5 


4.7 


12.3 


23.5 


10% 


D 


6.44 


1.00 


0.80 


7.4 


5.3 


15.0 


27.6 


4% 




0.75 


0.85 


1.38 




1.1 


19.0 


19.0 


17% 




2.47 


0.77 


0.78 




1.3 


16.0 


16.0 


30% 




4.35 


0.97 


0.80 




1.2 


12.1 


12.1 


58% 




6.35 


0.97 


0.78 




1.3 


15.2 


15.2 


35% 
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FIG. 2. Allowed regions at 95% C.L. for our effective two-neutrino fit to the solar neutrino 
data from Homestake, SAGE, GALLEX, and Super-Kamiokande. The four "finger" regions from 
Table are labeled A, B, C, and D. 
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FIG. 3. Allowed regions at 95% C.L. from the Super-Kamiokande atmospheric neutrino data 
for Sm^tm versus sin 2 202 when the overall atmospheric neutrino flux normalization a is allowed to 
vary, for sin$i = and sin#i = 0.2. 
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FIG. 4. Allowed regions at 95% C.L. from the Super-Kamiokande atmospheric neutrino data 
for Sm^ tm versus sin Q\ when sin 2 2#2 and the overall atmospheric neutrino flux normalization a 
are allowed to vary. The CHOOZ constraint is also shown. 
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FIG. 5. Allowed regions at 95% C.L. from the solar neutrino data for 5m1 un versus sin 2 29% 



when the B neutrino flux normalization (5 is allowed to vary, for (a) sin#i 
and (c) sin#i = 0.4 
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FIG. 6. Allowed regions at 95% C.L. from the solar neutrino data for 
when sin 2 2#3 and the 8 B neutrino flux normalization (3 are allowed to vary. 
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FIG. 7. Minimum \ 2 versus sin#i for the atmsopheric data set (dotted line), the solar data set 
(dashed), and the combined atmospheric and solar data sets (solid). All other relevant parameters 
are allowed to vary in these fits. 
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FIG. 8. Allowed region at 95% C.L. of the v e — > v T oscillation amplitude in atmospheric and 
long-baseline experiments versus 5m 2 atm , determined from a fit to the atmospheric neutrino data set 
when the CHOOZ constraint is included. The expected sensitivity of a long-baseline experiment 
from Fermilab to Soudan using muon storage rings is shown by the dashed line. 
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FIG. 9. Ratio of the electron energy spectrum to the SSM prediction for two different 
two-neutrino vacuum long-wavelength oscillation scenarios, compared to the Super-Kamiokande 
data. The shaded histograms show the results when oscillations are included. The two cases 
shown are (a) the best overall fit to all solar neutrino data, solution A from Table |TJ and (b) the 
best fit to just the Super-Kamiokande data, solution C from Table [□]. 
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FIG. 10. Electron energy spectrum in Super-Kamiokande for the Standard Solar Model with 
no oscillations (solid curve) and for the four vacuum neutrino oscillation solutions A (long-dashed) , 
B (short-dashed) , C (dotted) , and D (dash-dotted) listed in Table |TJ. Also shown is the current 
data from Super-Kamiokande [Q. 
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FIG. 11. Predicted event rate as a fraction of the SSM value in Super-Kamiokande versus time 
of the year for vacuum oscillation Solution A in Table || (dashed curve) and for no oscillations 
when the overall event rate is normalized to the observed value (solid). The other solutions in 
Table give similar results. Also shown is the current data from Super-Kamiokande |Hj. 
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FIG. 12. Predicted values for the seasonal asymmetry A\ denned in Eq. (^) versus electron 
energy in the SuperK experiment for the four vacuum neutrino oscillation solutions A (long-dashed 
curve), B (short-dashed), C (dotted), and D (dash-dotted) listed in Table ||. Also shown is the 
energy-independent value for no oscillations (solid). 
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FIG. 13. Typical predicted seasonal asymmetries A2 versus A± (denned in Eqs. ( |52| ) 
and (|53"|)) in the GALLEX and SAGE experiments for solutions from the four "fin- 
ger" regions of Fig. |: (a) sin 2 2<9 3 = 0.74, 5m 2 sun = (0.62 - 0.67) x lO" 10 eV 2 
and sin 2 2(9 3 = 0.91, Sm 2 sun = (0.73 - 0.77) x 10" 10 eV 2 (long-dashed curves), 
(b) sin 2 26*3 = 0.86, 5m 2 sun = (2.45 - 2.53) x lO" 10 eV 2 (short-dashed), (c) 
sin 2 26» 3 = 0.97, 5m 2 sun = (4.36 - 4.49) x lO" 10 eV 2 (dotted), and (d) sin 2 2<9 3 = 1.00, 
Sm 2 un = (6.38 — 6.51) x 10~ 10 eV 2 (dash-dotted). The predictions of best-fit solutions A, B, 
C, D from Table |TJ are indicated by the solid circles. In cases (c) and (d), 6m 2 un can vary over a 
wider range than shown; the additional values give orbits in A\— A% space that are slightly shifted 
from those shown here. The arrow on each curve is at the maximum value of 5m 2 un for each range. 
The no oscillation prediction (where the only variation in signal comes from the variation of the 
flux due to the elliptical orbit) is shown by a cross. 
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FIG. 14. Same as Fig. [T^ for the BOREXINO experiment. Cross sections have been calculated 
for final-state electron kinetic energies in the range 0.26 MeV < T e < 0.665 MeV. 
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